Non-Gaussianity from the hybrid potential 



O' 

(n: 



(N 






> 

en 
^, 

o 



X 



David Mulryne^'^, Stefano Orani^ and Arttu Rajantie^ 

^Astronomy Unit, School of Mathematical Sciences, 

Queen Mary, University of London, London, El 3NS. 

^Department of Physics, Imperial College, London, SW7 2AZ. 

(Dated: 21 July 2011) 

We study the hybrid inflationary potential in a regime where the defect field is light, and more than 
60 e-folds of accelerated expansion occur after the symmetry breaking transition. Using analytic and 
numerical techniques, we then identify parameter values within this regime for which the statistics 
of the primordial curvature perturbation are significantly non-Gaussian. Focusing on this range 
of parameters, we provide a specific example which leads to an observationally consistent power 
spectrum, and a level of non-Gaussianity within current WMAP bounds and in reach of the Planck 
satellite. An interesting feature of this example is that the initial conditions at horizon crossing 
appear quite natural. 



I. INTRODUCTION 

An inflationary epoch is the most promising model 
for the origin of large-scale structure in the universe |l| 
(for reviews see [2|), naturally accounting for current ob- 
servations which indicate that the universe's structure 
originated from near-scale invariant and almost Gaus- 
sian fluctuations [j, |j| . The small deviations from scale- 
invariance and Gaussianity are key probes of inflation, 
allowing competing models to be distinguished from one 
another by present and future observations. 

In canonical, single field inflation, the curvature per- 
turbation produced at horizon crossing is extremely close 
to Gaussian, and is conserved through the entirety of the 
subsequent evolution [5.-^]. On the other hand, when 
more than one field is light at horizon crossing, isocurva- 
ture modes are produced, which later can source an addi- 
tional contribution to the curvature perturbation. This 
can cause the curvature perturbation to become suffi- 
ciently non-Gaussian for this signal to be detected by 
future observations, or the model constrained by present 
ones. The current WMAP constraints on the /nl pa- 
rameter, which measures the deviation of the three-point 
correlation function from zero, is —10 < /nl < 74 at 95% 
confidence level, but /nl « 32 ± 21 at the 68% level [1]. 
Though premature, it is tempting to consider these num- 
bers as a 'tentative hint' of a non-Gaussian signal. The 
recently launched PLANCK satellite, will help resolve 
the issue. In the absence of a detection it is expected to 
give the constraint |/nl| < 5. 

With these issues in mind, considerable effort has been 
invested in developing methods to calculate how non- 
Gaussianity evolves on super-horizon scales during infla- 
tion, [9l4lll|. into identifying the features an inflationary 
potential must possess and the initial conditions required, 
to lead to a 'large' non-Gaussianity (/nl > 1) [13 - [l5| . 
Many scenarios have also been studied in which the gen- 
eration of non-Gaussianity occurs due to post inflationary 
effects, such as in curvaton ^6i], preheating [l7|, modu- 
lated (p)reheating ilSj], and models with an inhomoge- 
nious end of inflation [l^ [20| (see Ref. [2l| for connec- 
tions between these three scenarios). 



Of particular interest are effects which occur for infla- 
tionary potentials which are well motivated from a parti- 
cle physics perspective. The hybrid potential is one such 
example, and moreover is one of the most widely studied 
potentials in the literature. First introduced as an infla- 
tionary model by Linde ^^ , it is a two-field potential of 
the form 



^(0,x) = ^m2c/.2 



yc^'x' + ^ix'-vY- (1) 



As a quantum field theory, this is attractive because it 
is renormalisable and does not require unnatural cancel- 
lation of radiative corrections. In fact, potentials of this 
type are common in supersymmetric models [23|. 

In this paper we will investigate the hybrid potential 
^ using a different choice of parameters from the orig- 
inal proposal. We find that this leads to a viable cos- 
mological model, with perturbations that are compatible 
with current observations and level of non-Gaussianity 
that would be observable with the Planck satellite. 



II. ORIGINAL HYBRID INFLATION 

In the original version of the scenario [22| the (f> field 
plays the role of the inflaton. Initially, has a value 
larger than the 'critical' value, <j) > (/)crit = V^v/g, but 
still below the Planck scale < rupi = (SttG)"^/^. The 
latter requires Xv'^ < g^riA. The field then rolls down 
the potential towards smaller values until it reaches ^crit , 
at which point x becomes non-zero and inflation ends. 

We can generally consider a multi-field model with 
scalar fields </?" labelled by an integer a. In the case 
of the hybrid model, the fields are (fi = 4> and ip2 = X- 
Defining the slow-roll parameters. 
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the fields satisfy the 'slow-roll' equations of motion 
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The standard hybrid model assumes that the potential 
([1]) is dominated by the constant term during inflation, 
and this requires Aw* > m?m^y. The parameters are also 
chosen such that when </> reaches (/)crit, X rapidly grows 
and reaches the global minimum at % = ±u in less than 
an e-fold [2J|, thus ending inflation. This requires that 
the potential is steep enough in the x direction to avoid 
further inflation, or more precisely v < mp\. Also, 
must pass through (/)crit fast enough that x remains light 
for less than an e-fold. Utilising Eq. Q and © and 
linearising about 0crit such that A0 = </> — '/'crit, ^N — 
^e^^ijAc/), and m^ — 2g'^(j)critA(j). One finds that x is 
light for N w Au^/(m^?7ipj), and hence this combination 
of parameters must be less than unity. To summarise, 
therefore, the original scenario requires the parameter 
choices 



Xv^ < 
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V < TOpl, 

Xv < m m 
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It is important to realise, however, that these con- 
straints are not all necessarily required for a viable cos- 
mological model. They simply correspond to the specific 
physical scenario assumed in the original papers. De- 
spite having many appealing features, this original sce- 
nario suffers from two important issues. First, the phase 
transition generates a network of domains where differ- 
ent choices of vacuum have been made. The domains 
are separated by domain walls, which are incompatible 
with the observable universe. This issue can be over- 
come by using a multi-component x field, and possi- 
bly coupling it to a gauge field to eliminate Goldstone 
bosons, but this complicates the model significantly. Sec- 
ond, the scenario predicts a blue spectral index (n^ > 1). 
This is because for inflation driven solely by the (p field. 
Us = 1 — 6eJ + 2r?1 , (where * labels the time at which 
observational wavelengths became greater than the cos- 
mological horizon roughly 60 e-folds before the end of in- 
flation), and for the parameter choices discussed above, 
77! i, is positive and satisfies ryj , ^ e!. This is in conflict 
with WMAP and baryon acoustic oscillations measure- 
ments which give Ug = 0.968 ±0.012 with 68% confidence 
level il. 



III. HYBRID INFLATION WITH TWO LIGHT 
FIELDS 

Given the motivation for the form of potential (jlj , it 
is interesting to study its consequences outside of the pa- 
rameter ranges ([S]) which lead to the original scenario. 



For example, considerable interest has been devoted to 
parameter choices which make both fields light (see for 
example |25l,l27l|). and there have been a number of inves- 
tigations into the consequences of this situation, such as 
the resultant non-Gaussianity ^]. If the waterfall field 
is light at the symmetry breaking transition for more 
than an e-fold, however, the amplitude of density per- 
turbations on scales associated with this transition is in 
tension with primordial black hole constraints [23| . 

Recently, it has been noted by Clesse f2Sl that this con- 
clusion does not apply if the waterfall field is light but 
more than 60 e-folds of inflation occur after the phase 
transition. This is because any black holes produced are 
diluted by the expansion, and moreover, this also solves 
the defect problem, as defects are pushed onto scales 
larger than the present day observable universe. Clesse 
considered a parameter range identical to the original hy- 
brid scenario, except that he choose parameters such that 
(f) stays sufficiently close to (/)crit for x to be light for more 
than 60 e-folds, requiring \v^ 3> QQw^m^^. In this mini- 
mal adjustment to the original scenario, both fields still 
move sub-Planckian distances, but the velocity of the x 
field is much greater than the <j) field for the entire evolu- 
tion. This implies^ that n^ = 1 — Qe^ + '2-ri^x ^^^ because 
the model is still vacuum dominated and 7y^ is negative 
(because inflation occurs after the transition where x has 
a negative mass), a red spectral index results. The final 
issue with the original scenario is therefore resolved. 

Clesse's scenario does not lead to observable levels of 
non-Gaussianity. In this paper we show that there are 
other parameter choices which lead to perturbations that 
are compatible with observations but are significantly 
non-Gaussian. 



IV. THE 5N FORMALISM. 

To calculate the level of non-Gaussianity as well as the 
amplitude of the power spectrum and the spectral index, 
we will employ the 6N formalism [9, 2(^, |3^, implemented 
with numerical and analytic techniques. 

The 5N formalism, assumes the separate universe ap- 
proximation, in which points in the universe separated 
by more than a horizon size evolve like a 'separate uni- 
verse', the dynamics following from the local energy den- 
sity and pressure, and obeying the field equations of a 
Friedmann-Robertson- Walker (FRW) universe (see for 
example [3, [3l| ) . The 5N formalism then notes that the 
uniform density curvature perturbation between two such 
points, C, is given by the difference in e-folds of expansion 
between them, from some initial shared flat hypersurface 
labelled *, to some final shared hypersurface of uniform 



^ This is an approximation because in reality the dynamics of both 
fields contribute to the spectral index, as we shall see in the next 
section. 



density labelled /, 



C = SNi. 



(6) 



We use * to label the flat hypersurface, since in order to 
make observational predictions for an inflationary model, 
this surface must be defined at the time when obervable 
scales exited the cosmological horizon. Moreover, the 
final hypersurface / must be taken at some much later 
time when the dynamics are adiabatic and C is conserved. 
This can happen, for example, long after infiation ends, 
when the dynamics are dominated by a single fiuid. 

When slow-roll is a good approximation at horizon 
exit, ^ is completely determined by the field perturba- 
tions dip" at that time. 



C = <57V/(50. 



(7) 



The field perturbations are extremely close to Gaussian 
[33 |. and if their amplitude is sufficiently small, the statis- 
tics of the curvature perturbation can be determined in 
a simple manner by Taylor expanding 
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where here and from here on we employ the summation 
convention, and subsequently we will employ the notation 
in common use, Na = 9iV/((/'*, x*)/9(yc". 

One then finds that the amplitude of the power spec- 
trum is given by 



Al = No^N^ 
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which generalises the single field expression we employed 
earlier, and the amplitude of the reduced bispectrum by 
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V. ANALYTIC ESTIMATES 

In the end we will analyse the model using numerical 
methods, but to identify suitable parameter ranges, let us 
first use analytic approximations to model the behaviour 
of the model. One particular case in which analytic cal- 
culations are possible is when the potential is separable. 



non-separable term g^0^x^ is sub-dominant to the other 
terms. We will therefore use this approximation, and 
later check its validity using full numerical simulations. 
For separable potentials, there exist analytic expres- 
sions for the derivatives of N [l^, [SJl . The starting point 
for these formulae are the slow-roll equations of motion 
Eq. dH). It follows that 
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Therefore, for example. 
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The complexity of the calculation then reduces to deter- 
mining the partial derivative of the final field value on a 
constant energy density hypersurface with respect to the 
initial field value. However, when the dynamics become 
adiabatic during inflation, this derivative tends to zero. 
Therefore, one usually finds that^ N^ = ^x/('^piKxl*)- 
The second derivatives of N then follow by simple differ- 
entiation. 

One known mechanism which gives rise to a large non- 
Gaussianity in multiple field inflation is for one of the 
fields to be very close to a maximum of its self- interaction 
potential at the time of horizon crossing [IJ, [l^ , and for 
the magnitude of its tachyonic mass-squared to be much 
greater than the vacuum energy at this hilltop. In this 
region the x part of the potential can be approximated 
by Vi^ = Vo — 1/2to^x^i and the condition required is 
that rri^rn^^ ^ Vq [15]. Of course the field must also 
be light at horizon crossing, and hence rn^ < H^, but 
this condition can be simultaneously satisfied assuming 
the potential is not vacuum dominated, and that another 
field contributes significantly to the total energy density. 

Following Ref. [15[ we can understand this condi- 
tion by considering that when m^x^ < Vq, N^ « 
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i'ix*)j which becomes large in the limit that 
0. In this limit, the contribution to SN from 
field likely dominates over all other contribu- 



tions, and moreover iV^ 
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(5m2jm2)/(6yo). 

For the hybrid potential well after the symmetry break- 
ing transition, the shape of the potential is exactly of the 
form discussed above, with the non-separable term sub- 
dominant. In that case, Vq = Xv^/A, 
hence for this regime we expect 



Xv and 
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Of course, the hybrid potential Eq. ([T]) is not separable, 
but this should still be a reasonable approximation if the 



^ It is possible that this is not the case if Vx/V,x\f diverges as the 
derivative tends to zero (see [14|| for a full discussion). 



assuming that x* is sufficiently small. 

In principle, therefore, for v < rripi it is possible for 
this hybrid potential to give rise to a large positive 
non-Gaussianity. The additional constraint that the x 
field be light, however, introduces the restriction that 
l/2r7i^0^ ^ Vo = A'y*/4, and the part of the potential 
dominates the energy density. This implies we must move 
beyond the original parameter choices and consider large 
field values of in order that sufficient inflation occurs. 
This in turn implies that 0crit ^ ISmpi, where equality 
leads to roughly 60 e-folds after the transition. Hence 
4>* « 16mpi. In this case N^ = (/)*/(2mpj), and the con- 
dition that the x field sources the dominant contribution 
to (57V, lA'^^l > jiV^I, implies that v^ > 32x*mpi. We ar- 
rive at the curious situation in which the (/) field sources 
inflation, while the contribution due to the x field dom- 
inates C, and can be extremely non-Gaussian. This can 
be summarised by the following requirements on the pa- 
rameters 



Aw^ > g^mpi, 

m2,/)*2 > A^;^ 

V < mpi, 

0, sa 16mpi, 

X* ^ 0.03w^/mpi 



the numerator of Eq. (|Tl1) is still dominated by Xj we find 



(15) 



On the other hand, consideration of Eq. (ITU)) implies that 
the spectral index is generically red tilted and close to 
scale invariant, as will be confirmed by the full numerical 
simulations which follow. 

Guidance from analytic arguments, therefore, appears 
to lead to a large non-Gaussianity in the following sce- 
nario. Long before the last 60 e-folds of inflation, (j) be- 
gins its evolution with (j) > 0crit, and some value of x- 
Since the mass of x is initially positive it will evolve to- 
wards its minimum at x = as evolves towards 0ciit- 
Assuming the classical evolution of x reaches x = be- 
fore (j) reaches 0crit , then at this time the value of x will 
be dominated by its quantum diffusion during the pe- 
riod before the transition for which it was light. After 
the transition for the parameter choices discussed, more 
than 60 e-folds of inflation occur and x is still light when 
(j) = 16TOpi roughly 60 e-folds before the end of infiation. 
Its vev at this time will then likely still be close to x = 0. 
The X field will only roll significantly when x becomes 
heavy. If this occurs before cf) evolves to its minimum, 
it will roll during the slow-roll inflationary phase, oth- 
erwise will begin oscillating and begin to decay into 
radiation before x begins to roll. In the later case x will 
perhaps source an extremely short secondary inflationary 
phase as it rolls. It is clear, however, that Eq. P^ will 
remain at best a rough estimate. This is why we employ 
numerical simulations in our study. 

For completeness we note that when x* becomes larger 
and Eq. p^ becomes obsolete, it is still possible to find 
an analytic estimate for the asymptotic value of /nl- At 
X* ''^ v'^/rupi, we reach a regime where \N^\ > \N^\. If 
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indicating that /nl drops sharply. At even larger val- 
ues of x*7 we find N^j,^ > N^^, which implies /nl = 
5N^4,/{6Nl). Using N^ = <l>^l2ml^ and A^^ = 1/2™^^^ 
we obtain 
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Hence, we should find a smooth but clear transition from 
the significant non-Gaussianity predicted by Eq. ([T^ to 
practically Gaussian single-field behaviour ([T7| at around 
X* ~ 0.03w^/mpi. 

VI. INITIAL CONDITIONS 

The scenario discussed above requires a small initial 
value for the waterfall field, x* ^ v^/'Tipi, at </> = (/icrit- 
In this section we discuss the naturalness of this condi- 
tion. This is an instance of the 'measure problem' which 
is unresolved in general. However, we can apply some 
heuristic arguments. 

In Section|Vl we outlined a plausible sequence of events 
in which, long before the hybrid transition, both and x 
took large values, and the classical evolution of x tended 
to zero before </> reached (\>cnt- The initial values of (^ 
and X can be motivated by assuming a phase of eter- 
nal inflation during which the path the fields follow is 
dominated by quantum fiuctuations, rather than classi- 
cal rolling. The condition for such behaviour is that the 
distance moved in field space in a Hubble time, is less 
than a typical quantum fluctuation |34l | 



9 ^X^ > 



H' 
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(18) 



Using the slow roll equations of motion (U) this becomes 
a constraint on the fleld values, and deflnes a one di- 
mensional surface in field space at which the dynamics 
becomes predominantly classical. 

Beginning at this surface, therefore, one might won- 
der whether x will indeed reach x — (according to 
its classical rolling) before (f> reaches ^crit • The answer to 
this question will be extremely parameter dependent, and 
will likely also depend on the position on the surface from 
which the fields originate. It is easy, however, to build up 
a rough picture of the expected behaviour. Considering 
cases where both (j) and x are significantly displaced from 
zero initially, x will likely initially be the more massive 
field for the parameter choices we are focusing on (fT5|) . 
This is because x feels a mass of l/2Ax^ + g^(j)^ in com- 
parison with m? -\- g^x^, the mass of (/). Moreover, we 
expect X > g^ , unless v is much less than m^i. In the 
following we will assume this condition on A and g. The 





FIG. 1: The non-Gaussianity parameter /nl as a function of 
Xi and V, with other parameters given by Eq. (|2H) . Lighter 
shades indicate higher /nl- High values /nl ^ 100 are ruled 
out by WMAP, and small values /nl ;$ 1 are likely to be 
unobservable. 



FIG. 3: The ratio x*/Ht as a function of xt and «, with other 
parameters given by Eq. (|2H) . Lighter shades indicate higher 
X*/H,. The thicker contour hne corresponds to x*/H* — l. 
Values below this are unnatural because quantum fluctuations 
would generally give x* ^ ^* 




FIG. 2: The spectral index Us as a function of Xi and v, with 
other parameters given by Eq. (|2ip . Lighter shades indicate 
higher Us. All values obtained are compatible with WMAP 
observations. 
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FIG. 4: The dimensionless parameters A as a function of Xi 
and V, with other parameters given by Eq. H2ip . Lighter 
shades indicate higher A. The thicker contour line corresponds 
to A = IQ-i*. 



VII. NUMERICAL ANALYSIS 



mass of (j) would only be greater initially, therefore, if 



7,2 > 



\xi 



(where subscript e represents values on 



the eternal inflation surface). 

Assuming that x is indeed more massive initially, it 
will evolve rapidly towards x = Oi while (j) will remain 
nearly frozen. As x decreases, there will come a time 
when Xx^ < rrv^- At this point if the condition m^ > 
g^(j>^ is met then the x field will cease to be the more 
massive, and its evolution will slow significantly, and it 
will likely not reach x = 0. On the other hand, if this 
condition is not met we expect that x does reach zero. 
The necessary condition for x to reach zero, therefore, is 
that TO^ < g^0c, which ensures that x is always the more 
massive field. 

This condition is likely only very approximate, but 
from some limited numerical probing it appears to cap- 
ture the rough behaviour of the fields. In Section FVlIII we 
will look carefully at a particular parameter choice and 
probe the surface defined by Eq. (fT5|) . We find a com- 
plicated picture, however this rough analysis is a useful 
aid. 



Having identified the interesting parameter range, we 
now calculate the statistics of the perturbations numer- 
ically. Wc evolve the full non-slow-roll equations of mo- 
tion 
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and we have introduced a small decay term Ta to describe 
the decay of the scalar fields into radiation. As long as 
the decay rate is fairly slow, the produced perturbations 
are largely independent of it, so we do not attempt to 
use a realistic value but only check that our results do 
not depend on it. 

For simplicity we reduce the number of parameters 
which can be varied, by fixing 
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Xi/v < 1. 



(21) 



(22) 

Our parameter choice effectively means that the mag- 
nitude of iJ* and of T is proportional to vA, and so can 
be scaled simply by changing the value of A. On the 
other hand, since the derivatives of Ni , which enter for- 
mulae for observational quantities, are unaltered by such 
a scaling, changing the value of A leaves these derivatives 
unaltered. For each choice of Xi and v, we fix A such 
that the resulting amplitude A^ is in agreement with the 
observational requirement. The rescaling properties just 
discussed, however, mean that this can be done retro- 
spectively once the derivatives of N have been caculated 
(see below). 

As discussed in Section IIV[ we need to calculate the 
amount of expansion between a flat and an equal-energy 
density hypersurface, Nl {(f>^,x*), which we obtain by 
solving Eqs. (J19p for different initial values </>* and x*- 
The final hypersurface needs to be chosen in such a way 
that the scalar fields have decayed into radiation, so that 
the curvature perturbation is conserved. To find a suit- 
able final energy density or, equivalently, final Hubble 
rate, we first run the code from the initial condition 
4>i — 0crit and a suitable value of Xi- We follow the 
evolution until 99.9% of the energy density was in the 
radiation component, and use the corresponding Hubble 
rate Hf to define the final hypersurface. We also record 
the values 0* and x* i which the fields take 60 e-folds be- 
fore the end of infiation. The numerical system is then 
re-run from initial conditions which are slightly different 
to the recorded 4>^, and x* ■, and the number of e-folds to 
the same final value of Hf of the Hubble rate recorded, 
in such a way as to build up a finite difference approxi- 
mation for the derivatives of Ni {(j)^,, x*)- 

Once we have obtained Ni {(f>^,,x*), we calculate the 
spectral index Ug and /nl using Eqs. (fTO)) and (|TT|) . re- 
spectively. The results in Fig. [5] show that for all the pa- 
rameters we considered, the spectral index is contained 
in the region 0.9649 < Us < 0.9701. Therefore they are 
all compatible with the observations. 

On the other hand. Fig. [1] shows that a wide range of 
values for /nl are obtained, 7 x 10~^ < /nl ^ 5 x 10^. 
All parameters giving /nl ^ 100 are excluded by WMAP 
measurements [J|, constraining the parameter space. For 
Xi ^ 0.1f^/r?T,pi, values v < 0.07mpi are ruled out, and 
larger values of Xi are not constrained. 

Furthermore, Fig. [T] confirms the validity of the analyt- 
ical arguments and estimates in Section |V] . To a good 
approximation, we can identify x* ~ Xi- Our results 
show that as long as Xi is small enough, /nl is indepen- 
dent of it and decreases by one order of magnitude for a 
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FIG. 5: A'^ as a function of x* for the parameters shown in 
Eq. (|2H) . The result in Section [TVl assume that the quadratic 
Taylor expansion in Eq. ((5)1 is a good approximation. This 
is clearly not the case for the whole range of x*i but for a 
range that corresponds to the values present in one comoving 
volume corresponding to the currently observable universe it 
is. This is demonstrated by the inset, which shows the data 
for such a range with the linear term subtracted, together 
with a quadratic fit. 



two orders of magnitude increase in v/rripi, as predicted 
by Eq. pi]) . Towards the right, at Xi '^ O.lw^/rripi, /nl 
starts to fall in accordance with Eq. ([T6|) . and reaches 
a constant value /nl ~ 0.007 in good agreement with 
Eq. (HZD. 

As mentioned above in this section, with the 
parametrization given in Eq. (1211) , we can fix the parame- 
ters of the model retrospectively, constrained by the am- 
plitude of the power spectrum. As can be seen in Fig. 21 
10"^ < A < 10-25. 

An issue of importance is whether quantum fluctua- 
tions will affect which side of the potential the waterfall 
field is rolling along. This happens if the range of quan- 
tum fluctuations is greater than the distance of % to its 
hilltop. It is problematic because such a scenario leads 
to a configuration with domain walls, which we want to 
avoid. For a massless scalar field, the range is given by 
i/*/(27r). For a massive field, H^,/{2tt) is an upper bound. 
Thus, excluding parameters such that x*lii* ^ 1 settles 
the issue. As Fig. [3] shows, the excluded values are ap- 



proximately w < 5 X 10 ^rupi and Xi < v^ /fn-p\ 



VIII. A WORKING EXAMPLE 

Our results indicate that there is a fairly large window 
of parameters that are compatible with observations but 
lead to observable non-Gaussianity. We pick pick one 
representative example for more detailed study: 
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FIG. 6: Comparison of two initial conditions Xi ^nd Xi + ^Xi 
with parameters in Eq. (|21|) and Sxi ~ 1-4 x 10~'^77ipi. The 
plotted function is the difference in A*' as a function of en- 
ery density p, and corresponds to the curvature perturbation 
between points with the two different initial conditions. The 
vertical dashed line indicates when the waterfall field x starts 
rolling down the ridge. 



These parameters give 



/nl = 24.14, 
71, = 0.968. 



(24) 



Besides a value of the spectral index in the centre of 
the permitted observational range, these parameters give 
us a large, but observationally compatible, value of /nl- 
Indeed the value is close to the favoured WMAP value, 
/nl « 20 ± 10. 

In Fig. [5l we show N for a wide range of field val- 
ues X*- The expressions (fTTI) and (fTOj) for /nl and Ug 
assume quadratic Taylor expansion (|5]) of this function. 
It is clear from the figure that this assumption is not 
valid over the whole range shown. However, the actual 
range of Xi present in a comoving volume correspond- 
ing to the universe observable today is smaller, roughly 
6xi ~ y^N/2'KH^,, where iV « 60. The inset shows this 
range of x* , with the linear term subtracted. The remain- 
ing contribution is small and almost exactly quadratic, 
demonstrating that the Taylor expansion ([S]) is valid. 

It is interesting to confirm our expectation that the 
dominant contribution to 6N comes from the rolling of 
the waterfall field from its hilltop. This can be seen 
clearly in Fig jS] which plots the evolution of SN together 
with the point at which the waterfall field's evolution 
becomes significant. 

Finally, we ask whether the initial conditions we have 
assumed in this section are reasonable. Importantly, 
for the parameter values given in example (1231) . we find 
that the condition, given in Section IVH for mass struc- 
ture to be correct for x to evolve classically to x = 0, 
(/)^ > rn^/g'^, is met if 0c > 32mpi. This is necessarily 
satisfied when </> is above 0crit in this case. Of course, 
in practice we will need (j) to be significantly larger than 
this value in order that x ha-s time to evolve to zero. 
Numerically we find that (/)(, > ISOmpi is acceptable for 



FIG. 7: A plot of the boundary in field space on which quan- 
tum fluctuations become greater than the classical rolling of 
the fields. The three regions, described in the text, demarcate 
initial conditions on the surface for which x does not reach 
X = before (p reaches (jtait (far left), initial conditions which 
do (middle), and initial conditions for which eternal inflation 
does not end (far right). Because of the symmetry of the 
potential, only positive values of (j> and x a-re shown. 



X ^ v'^/iTipi at 0crit, which is satisfied on a significant 
proportion of the surface defined by Eq. (fTS]) . This can 
be seen by considering Fig. [71 which plots the surface on 
which eternal inflation ends for the parameter values at 
hand. For values of S> 150mpi the evolution is pre- 
dominantly in the x direction. Since the surface is com- 
plicated, an evolution which leaves the eternal regime 
may re-enter it, as can be seen from the figure. The 
dashed lines demarcate three regimes, trajectories origi- 
nating from the surface in the region on the far left do 
not evolve to x = before (j) reaches 0crit, while trajec- 
tories in the middle region do. Trajectories in the region 
on the far right originating from the upper eternal in- 
fiation boundary evolve classically towards x = Oi but 
reenter an eternal inflationary regime. The vast major- 
ity of trajectories which remain classical, however, evolve 
to X ^ v^/mpi when (j) ~ 4>cTit- 



IX. CONCLUSIONS 

In this paper we have demonstrated that inflation 
sourced by the hybrid potential Eq. ([T]) can give rise to 
an observationally significant non-Gaussianity while sat- 
isfying WMAP constraints on the power spectrum. The 
conditions required for this are summarised in Eq. (J15l) . 
This is an example of the hilltop mechanism for produc- 
ing non-Gaussianity first identified for axion models in 
Ref. lis]. 

In contrast with original hybrid infiation, where the 
waterfall field is heavy and cannot contribute signifi- 
cantly to the curvature perturbation on large scales }3a] , 
we take the waterfall field to be light. Our scenario also 
leads to more than 60 e-folds of accelerated expansion af- 
ter the waterfall transition. This avoids issues of primor- 
dial black holes and domain walls. While earlier work has 
considered the hybrid model with a light waterfall field 
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and the resulting non-Gaussianity [26|, l36| (and hybrid 
inflation with two hght inflaton fields [ij, [20]), we be- 
lieve that our work is the only study of non-Gaussianity 
in the case in which more than 60 e-folds of evolution 
occurs after the hybrid transition. 

Numerically we have only exhaustively studied the 
range of parameters discussed in Section Fyill but we have 
also probed several other parameter choices which lead to 
60 e-folds after the hybrid transition, including a number 
which are outside of the regime summarised by Eq. (|15p . 
Inside this regime analytic expectations are confirmed re- 
markably well by our simulations, which suggests that we 
should be able to trust them even outside it as long as 
the assumptions thay are based on are satisfied. Outside 
of this regime our non-exhaustive probing has not found 
any parameter choices which lead to a significant non- 
Gaussianity. It would be interesting, however, to probe 
an extended parameter range more carefully in order to 
determine whether the conditions Eq. ([T5|) are the only 
ones which lead to a signal. 

Finally we note that as discussed in Section I VII an in- 



teresting feature of the hybrid model is that it offers an 
explanation for the initial conditions at horizon crossing 
which lead to a large non-Gaussianity. While parameter 
choices are necessarily fine tuned, therefore, the initial 
conditions need not be, at least for some subset of the 
parameters which can produce a large non-Gaussianity. 
This is in contrast to other two field models which pro- 
duce a large non-Gaussianity due to infiationary dynam- 
ics [13 - [l4| . which require both carefully choosen model 
parameters, and offer no explanation for the origin of the 
finely tuned horizon crossing conditions. 
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